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Abstract 


The constrained minimization algorithm of Vanderplaats is applied 
to the problem of designing minimum drag faired bodies such as fuselages 
and nacelles. Body drag is computed by a variation of the Hess-Smith 
code. This variation includes a boundary layer computation. The encased 
payload provides arbitrary geometric constraints--specified a priori 
by the designer--below which the fairing cannot shrink. The optimiza- 
tion may include engine cooling air flows entering and exhausting 
through specific port locations on the body. 


Introduction 


Despite many protestations to the contrary, flight vehicle design 
has always been more art than science. It has, of course, been possible 
for some time to design simple shapes like projectiles according to 
formulae which can be shown to yield minimum wave drag for a given maxi- 
mum cross section and volume; but, for more complex shapes the process 
has always involved considerable cut and try as well as design experi- 
ence and aesthetics. 

However, this situation may soon change as the result of the con- 
fluence of three factors: (1) the development of practical algorithms 

to minimize objective functions which are subject to a variety of 
physical constraints; (2) the development of flow analysis techniques 
which can provide accurate results even for complex bodies at the end 
of which there is a small amount of flow separation; and (3) the arrival 
of computers with sufficient speed and memory capacity to run such 
codes efficiently and at reasonable cost. Before long it should be 
possible to routinely encase almost any size and shape payload in a 
minimum drag fairing specifically designed for that payload. It will 
also be possible to include in the optimization consideration items 
such as fairing weight, flow quantity into and out of the fairing, 
and the location on the body of the inlet and exhaust flow ports. 
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There are two very powerful economic incentives for following such 
an approach: (1) except for a small amount needed to increase the 

vehicle's energy, all of an aircraft's fuel is burned overcoming drag 
and (2) a cut and try approach to aircraft drag minimization is ex- 
tremely costly in time, in engineering and technician labor, and even 
to some extent in wasted materials. 

It has long been recognized that the drag increment due to the in- 
stallation of an air-cooled piston engine is related in some fashion to 
its cross sectional area. It has, therefore, been the practice to cowl 
it as tightly as possible while leaving adequate room for accessories, 
ducting, and cooling air flow. Within these constraints it is usually 
left to the designer's eye to determine the exact cowl shape. Only if 
the resulting aircraft's performance is far below predictions would one 
normally expect to examine the adequacy of the nacelle or forward fuse- 
lage design by way of test program perhaps combined with an extensive 
analytical investigation. 

It is also known that the drag associated with the flow of cooling 
air through nacelle inlets, engine fins and baffles, and nacelle ex- 
hausts can represent as much as 25% of the total vehicle drag.^ This 
figure, even without considering the drag increment due to an improperly 
shaped cowl, is sufficiently large to justify serious attempts at de- 
veloping a rigorous design technique guaranteed to produce a minimum- 
drag, nacelle-flow combination. 

An analytical optimization technique is desired which can indicate 
how the body shape should be modified so as to yield a demonstrable. . 
minimum drag. Obviously, this must be coupled to a sufficiently accu- 
rate fluid dynamic model of the vehicle so that the measured drag re- 
sulting from these shape changes will in fact be very close to the 
drag calculated. This is a report on a NASA-supported research project 
at North Carolina State University which has this objective as its 
goal . 


The Optimizer 


The function to be minimized is the steady-state drag of a fairing 
about a specific "payload." Payload is here taken to mean such things 
as an engine, a passenger compartment, etc. The payload provides 
specific geometric constraints limiting the minimum fairing size. 

The payload is not necessarily regular or symmetric about any 
axis. A variety of optimizers could be used. It has, however, been 
the experience of many workers in the optimization field that the best 
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method to use in a specific situation often depends upon the problem at 
hand and cannot readily be foretold. Those who have worked in the 
field for some time are also aware of the difficulty of developing new 
codes--even those following established procedures --from scratch. The 
use of generally available optimization codes whenever possible is 
therefore highly desirable. One of the most versatile of these is that 
developed over the past ten years or so by Garrett N. Vanderplaats^ of 
the NASA Ames Laboratory and the Naval Postgraduate School. Further, 
it has been applied successfully to a nunier of related optimization 
problems (airfoil design). 3 This code employs the method of feasible 
di recti ons. 

The Vanderplaats procedure has been formulated as a numerical 
optimization problem which may be stated as follows: 

minimize OBJ = F(X) 
subject to G^(X) < p (i = l,m) 

xf < X. < (i = l.n) 

In the present case, the objective function is the body drag. The vec- 
tor X contains n design variables. For axisymmetric bodies we have 
chosen n = 5. For more general bodies with a plane of symmetry repre- 
sented by a set of 14 panels over the half body, we take n = 15 * 5 75 

The vector J contains the amplitudes of the shape functions (0 < x _< 1) 


(Here x refers to the non- 
dimensional distance from 
the body nose.) 


which are employed to change the shape of the body locally along a line 
reaching from the nose to the tail.* For the axisymmetric body, every 
ci rcumferential panel is given the same modification when one shape 


F , yO.ZS * 


20 X 


F. = 


F. = 


F, = 


fr. = 


sin (. x°-^ 31^3 
sin (. 

. , 1.357.3 

sin (it X ) 

. / 3.106,3 

sin (it X ) 


*The exponents on x determine the axial station at which that function 
produces its maximum influence on the body shape. These exponent values 
may be problem dependent. 
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function amplitude is varied. For the general body, the surface shape 
modifications due to a change in one shape function amplitude are 
applied along a line of panel boundaries reaching from nose to tail. 

The shape function amplitudes in J are to be changed in such a way that 
drag is minimized satisfying the constraints. Gj(X) defines the con- 
straints_which thejdesigner wishes to impose on the optimization prob- 
lem. F(X) and G^(X) may be either implicit or explicit functions of 
th^ design variables X b^t must be continuous. In the present problem, 
F(X) is implicit but G^-(X) is explicit, representing the m major physi- 
cal dimensions of the inscribed payload. 

— ‘Q — ^ 

Variables X^ and X^- define the lower and upper bounds, respe^c- 
tive^ly, on the design variables and are the limits over which F(X) and 
G<(X) are defined. For the pro^em being considered is unbounded. 

If the inequality condition > 0 is violated for any constraint, 

the_constraint is said to be violated. If the equality condition G^ 
G^-(X) = 0 is met, the cojistraint is called active and if the strict 
inequality condition G^-(X) < 0 is met, the constraint is inactive. 
Computationally an exact zero can seldom be aclrieved; thus a constraint 
is called active if its value is within a specified tolerance of zero. 

The n-dimensional space spanned by the design variables X is re- 
ferred to as the design space. Any design which satisfies the in- 
equality equations is referred to as a feasible design. If the design 
violates one or more of these inequalities, it is said to be infeasible. 
The minimum feasible design is said to be optimal. 

The optimization program begins with an initial X vector which is 
input to the program and may or may not define a feasible design. (In 
the present case, we always begin with a circumscribing body larger 
than the minimum drag body and therefore always feasible design.) The 
optimization process then proceeds iteratively by the following recur- 
sive relationship: 

j?’*' = jTi + c* 

where q is the iteration nuirber, vector S is the direction of search 
in the n-dimension design space, and a* is a scalar which defines dis- 
tance of travel and direction S. The notation a* for the move para- 
meter is used here for consistency with mathematical programming 
nomenclature and should not be confused with the body angle of attack 
which, for reasons indicated below, is always zero. 

The optimization process then proceeds in two steps. The first 
is the determination of a direction S which will reduce the objective 
function without violating the constraints. The second is the deter- 
mination of a scalar a* so that the objective function is minimized in 
this direction, a new constraint is encountered, or a currently active 
constraint is encountered again. For exanple, suppose an initial 
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design is given for which no constraints are active or violated. The 
program then begins by perturbing each of the X variables to determine 
its effect on the objective function (body drag, in this case). The 
gradient of Cq is calculated by finite difference using a single forward 
step and the gradient vector is constructed as 


VOBJ 


Because no constraints are active 
greatest improvement in the objective function is obtained by mgving_ 
in the negative gradient or steepest descent direction so that S = -VC^ 

Having determined S, the scalar a* must now be determined so that 
either the objective function is minimized in this direction or some 
constraint boundary i£ encountered. That is, one-dimensional search 
is done in direction S to determine the appropriate value for a* so 
that an improved design is obtained. No further improvement can be_ 
achieved in this direction; it is now necessary to determine a new 5 
vector which will improve the design. 

The second design iteration begins by again perturbing the design 
variables to obtain the gradient of the objective function, Cq. Now, 
instead of moving in a steepest descent direction, the S vector is cal- 
culated from the following relationship 



or violated, it is obvious that the 


2 

^ . .tobjo + sq-1 

I 

This equation defines the conjugate direction as developed by Fletcher 
and Reeves. The advantage that this definition of the S vector has over 
the steepest descent direction is that if the objective is a quadratic 
function of the design variables, covergence to the optimum can be 
guaranteed to occur in n iterations or less. Although most problems 
of practical interest are not quadratic functions of the design varia- 
bles, they may still be approximated as quadratic in the region of 
the solution. In other words, if the first three terms of a Taylor 
series expansion of the objective function form a reasonable 
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approximation to that function, then the foregoing equation can be ex- 
pected to provide more rapid convergence than a steepest descent 
search, since steepest descent uses only the first two terms of the 
Taylor series expansion. Having determined the new direction S, one 
searches in that direction until another constraint is encountered, 
which ends the second design iteration. 

The design variables are again perturbed to obtain the gradient of 
the objective function, C^. At the same time, the gradient of the ac- 
tive constraint is obtained. Now a search direction must be found 
which will reduce the objective function without violating the active 
constraint. Such a direction can be found by solving the following 
subproblem which is a linear programming problem with a single quad- 
ratic constraint. 

Find S to maximize 6 ^ 

Subject to: 

VOBJ (X ) • S + 6 < 0 

• S + B <0 (j = 1,NAC) 

S • S < 1 

MAC is the nunter of active constraints. The details for solving this 
problem are given by Zoutendijk^ and by Vanderplaats and Moses. ^ Note 
that if the first equation is satisfied and B is positive, the result- 
ing direction will reduce the objective function and is defined as a 
"usable" direction. If the second equation is satisfied and B is posi- 
tive, S is called a feasible direction because, for a small move in 
this direction, no constraints will be violated. The prespecified 
parameter 0j is referred to as a pushoff factor for the jth constraint 
and has the effect of pushing and design away from the active con- 
straint. The value of 0j must be zero or positive to maintain a feasi- 
ble design. If 8j were zero, the resulting direction would be precisely 
tangent to the active constraint. On the other hand, a very large 9j 
would push the design away from the active constraint and nearly tan- 
gent to a line of constant objective function. A value of 0j = 1 will 
yield a direction which approximately bisects the angle between con- 
stant objective function and the constraint. If the maximum value of 
B obtainable from the above equations is zero, then no direction exists 
which will both reduce the objective function and s-atisfy the constraint. 
The current design is optimal or at least a local minimum. 



ORIGtNAl PAGE IS 
OF POOR QUALITY 

7 


As the number of design variables increase and the nonlinearity of 
the objective function increases, it will be necessary in practice to 
limit the size of a* to prevent the body fi’om becoming so distorted that 
that fluid dynamic analysis is not possible. Increasing the nunter of 
design variables and the nonlinearity of the objective function also 
increases the accuracy requirements on the fluid dynamic calculations 
since each design variable now contributes a smaller amount to the 
total gradient value and the contributions from different variables 
may range from nearly insignificant to very significant, a situation 
that may change from iteration to iteration as the body shape changes. 
Even without these accuracy requirements and move limitations, the 
computational time required increases approximately as the number of 
design variables squared. Thus, while a two-dimensional airfoil 
optimization problem involving 5 design variables may require perhaps 
two minutes per iteration and no more than 20 iterations to reach an 
optimum that also satisfies the constraints, the minimum drag body 
problem involving 75 design variables requires on the order of an hour 
per iteration and--depending upon the extent to which the optimization 
is carried out, e.g., whether the optimum is actually reached or merely 
closely approached- -many more than 20 iterations. 

The use of such mathematical techniques in design then becomes a 
question of economic trade-offs. The speed of digital processing is 
increasing and its relative cost is decreasing. Manual design labor, 
operation of test facilities for evaluating the effectiveness of body 
recontouring, and the energy cost of operating a body with more than 
optimal drag are all increasing. Thus it seems likely that if a 
mathematical analysis and optimization technique can be shown to yield 
accurate results, it will replace the traditional design and develop- 
ment technique in the near future. 


Flow Analysis Model 


In order to be useful in real situations, the flow model upon 
which the optimizer operates must be capable of accurately represent- 
ing the flow over rather complex bodies. The authors had been involved 
in the development of such a code during the previous six years so it 
was natural for them to consider that code as the basis of the flow 
model for the present work. As is usually the case in such matters the 
code has a long pedigree. It began in 1962 as the work of John Hess 
and A. M. 0. Smith^ of the Douglas Aircraft Company in Long Beach, 
California. They chose to represent the potential flow over a body 
with a plane of symmetry through the following device: the surface of 

the body is represented by many plane quadrilaterals. (In the limit 
^is representation becomes exact.) At the center of each panel is a 
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fluid source of undetemrined strength. (One could also use a distributed 
source of constant strength over a panel.) The source strengths are 
determined by requiring that the flow induced by all sources interacting 
with the oncoming stream be parallel with the surface at one point on 
each panel.* 

Charles W. Dawson and Janet Dean^ of the Naval Ship Research and 
Development Center added to the Hess-Smith code the capability for cal- 
culating on-body and off-body streamlines in the early 1970's. At North 
Carolina State University we used this capability to do two-dimensional 
momentum integral boundary layer calculations along streamlines. 8 As a 
result, we were able to determine the skin friction at the center of 
each panel and the boundary layer displacement thickness in the aft 
regions of the body. Based on rather detailed surface pressure measure- 
ments on the Airship Akron which we found in the literature, we chose 
to assume that the flow separated from the body at a point two sets of 
panels forward from the rear of the body with an angle relative to the 
surface proportional to d6*/dx. We then located the downstream stagna- 
tion point in such a way that the departing streamline would gradually . 
fair to it, about one body diameter aft of the physical body. For cal- 
culating the pressure distribution on the body we assumed that the phys- 
ical stem end of the body was replaced by this faired "wake body" and 
the pressures calculated for the first two sets of panels on the wake 
body were those applied to the last two sets of panels on the physical 
body. As a result, the real body evidences both a skin friction drag 
and a pressure drag. 

It will be recognized that this device for obtaining body drag 
within the framework of a potential flow calculation rests on a nunfcer 
of rather arbitrary assumptions: 

1. That the boundary layer flow can be analyzed as being 
two-dimensional along streamlines. This is approxi- 
mately correct only when we are dealing with a rela- 
tively slender body at zero angle of attack. The 
method obviously fails in situations with large cir- 
cumferential pressure gradients. 


*Inlets and exhausts are represented by relaxing the requirement that 
the flow velocity normal to every panel surface be zero. This require- 
ment is replaced by one specifying the direction and magnitude of the 
flow velocity normal to the surfaces of these panels which form inlets 
or exhausts the zero normal velocity requirement is retained for all 
other panels. 



ORiGfNAL PAGE IS 
OF POOR QUALITY 


2 . That a standard motnentum-integral technique may be 
used. Momentum integral techniques have been shown 
to be reasonably accurate for constant velocity and 
accelerating flows. They do not give an accurate 
indication of the point of flow separation and yield 
too large a boundary layer in regions of decelerating 
flow--hence, the necessity of specifying the separa- 
tion point (point for attachment of the wake body) 
rather arbitrarily. For this reason also, they give 
poor results for bodies with concave regions. (The 
potential flow method also has trouble with such 
regions.) Finite difference calculation of the 
boundary layer characteristics were rejected because 
of the very large computing time penalty involved 
and the present uncertain state of general three- 
dimensional boundary layer theory. 

3. That the point of attachment and shape of the wake , 

body can be specified a priori based on a limited '' 

correlation of body shape and drag. 

Despite these arbitrary assunptions, calculated drag values were 
acceptably close to experimental values (less than 5% error) for those 
fuselages and nacelles for which we could find experimental results. 

The version of the NCSU BODY code used for the present work does 
not include the panel normal velocity specification provision developed 
previously 5 which permits the representation of flow inlets and exhaust 
by specifying the entering and leaving flow velocities based on the ex- 
ternal pressures and the internal flew restrictions. It also does not 
include a provision, supplied in an earlier version'^, for imposing the 
axial component of propeller slipstreams on the flow over the body. 

These provisions were omitted for the present in order to hold down 
the cost of runs during code development. One other cost-saving change 
of a rather fundamental nature was to eliminate most of the disk read- 
writes present in previous versions. This step resulted in a 40% im- 
provement in execution time and a 30% reduction in cost. It also ren- 
ders the code too large to execute on most computers without virtual 
memory capabilities. 

Note, too, that the version of the potential flow code used here 
does not include some of the later provisions added by Hess and Smith 
such as the use of linear source distributions on curved panels. These 
provisions improve the accuracy of the code significantly for a modest 
increase in computer time. 
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Mating the Optimizer to the Flow Analyzer 


In order to determine how the body shape must be altered to reduce 
the drag, the optimizer must evaluate a gradient (A drag/A shape) nu- 
merically. This gradient is a vector which contains an eleirent repre- 
senting the contribution from each variable specifying the body surface 
If one were to use the comer coordinates of each panel as design 
variables one would have to go through the body flow analysis program 
N times, where N is the nunber of panels used to represent a half-body. 
When N is 600 or more, as it needs to be to adequately represent even 
moderately complex bodies, the conputer time required becomes imprac- 
tical even in this day of high speed computers. For this reason the 
use of shaping functions or reduced coordinates appears very attractive 
These alter the body shape with fewer design variables. Each function 
alters a small region of the body surface rather than just a group of 
four panels. Some flexibility is sacrificed but by judicious choice of 
the shaping functions virtually any modification felicitous for drag 
reduction can be implemented. For our work, we choose to employ the 
same shape functions Haney and Johnson^ 1 used for their airfoil opti- 
mization work. Using this approach, we employ a group of five shape 
functions to represent the entire length of the body at each circum- 
ferential panel boundary. In other words, if one has a body repre- 
sented by N axial panel boundaries and M circumferential panel bound- 
aries for a total of (N-1)*(M-1) panels one can, through the use of 
these reduced coordinates (shape function amplitudes), reduce the num- 
ber of design variables from N*M to 5*M. 

This, then, was the plan of attack. As with many other plans, 
difficulties not anticipated in the conception phase arose during 
execution. The computation of the gradient involves taking differ- 
ences between the drag of the unaltered body and that of the body with 
one design variable altered. Since the difference between the two 
values usually occurs in the third and fourth significant digits, an 
error of 5% in the drag computation is obviously not acceptable. This 
problem is evidenced by a rapidly varying sign on the gradient from 
iteration to iteration. Modifying the flow analyzer to execute in 
double precision seemed to resolve that problem. 

It was noted earlier that the wake body was attached to the physi- 
cal body two set of panels forward of the end of the physical body. 
While this arbitrary attachment point proved satisfactory for rela- 
tively streamlined bodies, the optimization scheme tended to squeeze 
the fore and aft ends of the physical body, leading to very steep 
pressure gradients toward the rear of the body. The computed drag 
would certainly decline as a result of such squeezing but in reality 
the flow cannot negotiate the large pressure rises indicated by the 
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computations in the vicinity of the rear end of the body. The obvious 
remedy is to move the wake body attachment point upstream to the point 
where the flow would normally separate under these circumstances and 
extend the downstream tip of the wake body (the stagnation point re- 
quired by potential theory) further into the physical wake so that 
the rate of pressure rise with streamwise distance on the physical 
body is more realistic. Determination of when it is time to perform 
this action can be made in a number of ways. As discussed in the next 
section, the approach used in formulating the present code is to test 
df/dx all along the aft portion of the body-wake combination. If this 
value ever exceeds a trigger value, it is assumed that separation would 
have taken place. The wake attachment point is then moved upstream 
and the stagnation point downstream and the pressures are recomputed. 

The resulting wake bodies tend to become very long, shallow-angle cones, 
particularly if the attachment point had to be moved more than once to 
obtain an acceptable pressure distribution. 

The existence of a large, gradually-tapered wake body can easily 
illuminate the existence of another problem. An integration of the 
pressures over the body cum wake body should yield, according to 
D'Alembert's principal, zero drag. However, the pressures are known 
only at the centroids of N flat panels representing the body surface. 

It is natural that one should carry out the force determination by 
assuming these pressures to exist over the entire panel and multiply 
the product of pressure and panel area by the cosine of the angle be- 
tween the panel normal and the streanviise direction and then summing 
over all panels. For a long, thin wake body the pressure at the last 
panel's centroid, for example, is quite different from that at the 
downstream stagnation point. For an elliptical body represented by 
400 panels errors as high as 35% can arise from this source. It is 
not practical to increase the number of panels used in the computation 
significantly since storage increases directly with N and the compu- 
tational time per iteration increases as n 2 while the error is only 
one-half for twice as many panels. It is, therefore, necessary to 
attempt to fair the pressure distribution smoothly over the body and 
to integrate the variable pressure over each panel. This is a rather 
formidable task which seems to be best acconplished by transforming 
the panel to a square and employing a standard bi-cubic spline routine* 
to accomplish the pressure fairing and force integration. 


*From the IMSL package. 
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Results and Pi s cussion 


Because of the shorter run times required, much of the development 
work with the code used the axisymmetric body option. The constraint 
body was a 3-to-l ellipsoid and the circumscribing body a 2-to-l ellip- 
soid. Figure 1 shows the ci rcumscribing body. Figure 2 shows the body 
with the wake body attached over the last two sets of panels on the 
physical body. The drag coefficient which this arrangement yields is 
indicated beside the body.* The value reported in Rouse 12 for a 2-to-l 
ellipsoid with turbulent boundary layer using the present normalization 
is 0.0127-0.01459. The body and wake configuration which the optimizer 
yields if one does not permit the separation point to move is shown in 
Figure 3. Figure 4 shows the wake and body combination obtained by 
allowing the separation point to move upstream in response to an exces- 
sive pressure rise near the aft end. This figure illustrates particu- 
larly well why the pressure at the center of the last set of panels — 
if assumed to exist over the entire panel— will integrate to a non- 
zero drag and why it is necessary to employ a variable pressure-area 
relationship, at least in this region, for accurate drag determination. 

Since the wake is the crucial factor in determining the drag (the 
skin friction is only about 20% of the total drag of the original 
circumscribing body and does not change significantly with small varia- 
tions in body configuration) the wake will also be the critical factor 
in determining whether and how fast the body diameter decreases to 
meet the constraint body. If the optimizer can shrink the wake alone 
it has little incentive to reduce the diameter of the circumscribing 
body. Thus, it is important that the gradient of the objective func- 
tion be capable of including the effects of moving the separation point 
upstream and that a realistic criterion for the initiation of this 
movement be included in the routine. 

At present, the routine calculates the change in pressure from 
panel centroid to panel centroid and divides this by the axial distance 
between centroids. Whenever this pressure gradient exceeds a present 
value (selected arbitrarily at this time by the user) it restarts the 


*Note that one can make this number almost anything one desires by 
altering the location of the downstream stagnation point and the way the 
the "wake" is faired to the body. A downstream stagnation point lo- 
cated one body diameter aft of the stern coirfcined with a wake having a 
rapid initial taper (so that the wake looks rather like a golf tee) 
will yield a low Cq value. On the other hand, if the downstream stag- 
nation point is located, say, 10 body diameters downstream and the wake 
has almost no taper, Cg will be rather large. A 3-to-l ellipsoid with 
a cylindrical wake attached at the maximum thickness station would 
yield a Cg of about 0.08 with the present method of normalization. 
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program with the wake attached two axial stations upstream. The pres- 
sure gradients are recomputed and if they again violate the preset 
limit, the wake attachment point is again moved upstream two more 
axial stations. The process is repeated until the pressure gradient on 
the body plus wake is everywhere below the preset value or the wake is 
attached at the axial station corresponding to the maximum body dia- 
meter or at the first station downstream of the maximum body diameter. 
Because the pressures on the aft end of the physical body ar^ then 
quite low, the optimizer wi 11 shrink the body rapidly as shown in 
Figure 5. 

If the "pushoff factor" is set too small, the body produced by the 
optimizer will violate the constraints. See Figure 6. On the other 
hand, making the constraints too elastic (pushoff factor too large) re- 
sults in the generation of a large disc-like body with a rapidly 
tapering wake once the constraints are encountered. This type of wake 
occurs because of the necessity of specifying the downstream extent of 
the wake a priori in terms of the original body diameter. For this 
reason the wake contours become very unrealistic if the diameter of the 
body produced by the optimizer is greatly increased. Figures 7 and 8 
depict this rather forcefully. 

In operating the program one must also learn to adjust the rela- 
tive strengths of the shaping functions and the axial station of maxi- 
mum application for each function. Figures 9 and 10 show the result of 
insufficient strength in the rear most shaping function. The forebody 
is well contoured and the overall body plus wake is quite streamlined. 
The physical body, however, is rather blunt toward the stern end. The 
maximum diameter is about 12% larger than necessary and the drag is 
about 50% larger than it could be. By increasing the relative strength 
of the shaping function positioned with its maximum effect at the rear 
end of the physical body, the body shown in Figures 11 and 12 is ob- 
tained. The forebody shape is about the same as before but the maximum 
diameter has been reduced to 1% larger than the maximum constraint body 
diameter and the aft end is now tririmer. While the body may still 
appear to be somewhat blunt, the conputed pressure drag coefficient of 
0.00031 compared with an overall Cq of 0.00254* indicates that any 


*Drag coefficient is based on an arbitrary reference area which is used 
for all computations. The changes in Cq therefore indicate the changes 
in actual drag resulting from body shape modifications. Ref. 12 gives 
0.0010 as the Cq for an airship hull when normalized in this fashion. 
This value represents the limit to which the present calculation should 
proceed if the circumscribing body is positioned such that a smoothly 
tapering tail of minimum length can be produced. For comparison pur- 
poses, note that the surface area. of a 6-to-l prolate spheroid with 
major axis equal to 12 is 59.94. (The final body has a calculated 
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additional body shape modifications can at best result in a 10% reduc- 
tion in drag** and that these modifications must be confined the last 
two sets of panels. To effect this change, the relative importance of 
the shape function which has its maximum influence in this region must 
be increased significantly. In retrospect, it appears that positioning 
the constraint body further forward in the circumscribing body would 
have provided additional panel sets for the optimizer to operate upon 
and would have yielded a more streamlined afterbody shape. 

As the program now exists it can modify the circumscribing body 
only in the Y and Z directions. The length in the streamwise direction 
is fixed. Thus, in using the program one should take care to locate 
the upstream and downstream ends of the circumscribing body at those 
locations which experience indicates will yield a reasonable final 
shape. The program will then modify the Y and Z dimensions to produce 
the minimum drag possible for these locations consistent with con- 
straint body dimensions. 

It will be recognized that in this program we have attempted to 
model an inherently viscous flow phenomenon through the device of rep- 
representing the separated viscous wake behind the body by a pseudo- 
surface. This then permits one to obtain results agreeing reasonably 
well with reality from an inviscid flow analysis. The reason one must 
resort to such tricks is that the analysis of a fully viscous flow 
about a body is a problem perhaps two orders of magnitude more diffi- 
cult than that of an inviscid flow about the same body. The strategem 
is not without pitfalls, however. Additional information of the type 
which would naturally be developed during the solution of the viscous 
flow problem must be supplied a priori to locate the wake attachment 
point and the downstream "stagnation" point. Since this kind of infor- 
mation is developed but slowly through the experimental correlations. 


surface area of 58.37.) A 3-to-l prolate spheroid with major axis 
equal to 6 has a surface area equal to 30.89. Replacing the aft half 
of a 3-to-l ellipsoid by a right circular cone 9 units long yields a 
body with a surface area of 45.25. Assuming that the frictional drag 
is directly proportional to the surface area and that the remaining 
pressure drag can be eliminated through the use of a conical tail, 
the limiting Cq for the procedure becomes 0.00167. A further reduc- 
tion would be possible if the tail could be shortened without induc- 
ing a pressure drag. 

**The skin friction cannot be significantly reduced because the over- 
all body length remains fixed and the minimum diameter is limited by 
a constraint body. 
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one must proceed cautiously. Nevertheless, the evidence to date indi- 
cates that astute modeling combined with optimal control techniques 
can result in a superior, cost-effective vehicle design technique. 

Organizations wishing to experiment with the code in its present 
form can obtain copies through the contract monitor. Because computer 
core is much less expensive at NCSU U\an frequent disc read-writes, 
the program is configured to take advantage of all the core that could 
be conveniently assigned to one job (about 6,000 K-Bytes on an IBM 3081). 
This allocation permits one to provide about 196 panels on the physical 
body plus 28 for the wake body. Execution time is about 8 minutes per 
iteration. In principle, when larger computer memories and higher speed 
processors become available, bodies of greater complexity can be 
treated merely by increasing the array sizes and the appropriate DO 
1 OOP 1 i mi ts . 

The authors will endeavor to answer the questions of potential . 
users of the code. A detailed user's manual is being prepared and 
will be available separately. This will include sample inputs, varia- 
ble identification, and sample outputs. 


Recommendation for Future Work 


A computer code of this magnitude is necessarily evolutionary. 
Some potential additions which the authors feel should be investigated 
for possible inclusion are: 

1. The panel normal velocity specification provision to per- 
mit the treatment of cooling flows. This is a relatively 
simple step, having been included in previous versions of 
the flow analysis code. It has been intended that this be 
done during the current grant period. Data sets for an 
actual light-twin engine nacelle and a possible redesign 
were prepared for use when this feature became available. 
Unexpected delays, principally those due to providing 
variable wake-body sizes and attachments , used the time 
allotted for this purpose, however. 

2. The provision for treating propeller slipstreams. The 
previously-used package which computes only the axial 
flow component can be added relatively simply but the 
authors feel that replacing this with a computation em- 
ploying trailing helical vortices, which can include the 
effect of wings protruding large distances from the 
body, should be strongly considered. 
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3. More rigorous separation criteria based on the velocity 
distributions in the local boundary layers for both 
laminar and turbulent flow. This would replace the 
rather arbitrary separation criterion now used which is 
based solely on the external pressure gradient. 

4. A provision for accommodating small angles of attack. 

The original XYZ potential flow program included a pro- 
vision for adding the effects of onset flows in the Y 
and Z directions. The streamline confutation can also 
accommodate these cross flows. If one assumes that the 
boundary layer is thin enough that all of it moves in 
the same direction as the external flow, then the exist- 
ing boundary layer computation method, modified in some 
fashion to account for the more rapid spreading of the 
streamlines associated with cross flows when compared 
with the axial flow over thin bodies and the more 
rapid curvature changes in the cross flow direction may 
possibly yield results which are adequate for design 
purposes if the velocity distribution in the neighbor- 
hood of the separation point can, at the same time, be 
determined more precisely. A complete finite-difference 
computation is not desired because of the computational 
time required. However, it may prove desirable to begin 
a two-dimensional finite difference computation at the 
point along each streamline where the pressure gradient 
becomes adverse. 

5. A provision for including an axial scaling of the cir- 
cumscribing body within the optimization procedure. 

Possibly this could take the form of two additional de- 
sign variables, one applied to front half of the body 
and the other to the rear half. 

The authors also feel that the user community would profit from 
additional exercise of the code. It would be possible thereby to ac- 
quire a better understanding of the effect of input parameter changes, 
location of the constraint body within the circumscribing body, and 
size relationship between the constraint body and the ci rcumscribing 
body. They would be pleased to work with organizations endeavoring to 
undertake this activity. 
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Figure 1. Circumscribing Ellipsoid 



Figure 2. Ellipsoid with Wake Body 
(Enlarged) 
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Figure 3. Ellipsoid with Coatracted 5ow_and Stern 



Figure 4. Ellipsoid with Extended Wake Body 
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Figure 5. Shrunk Body 



Figure 6. Body with Violated Constraints 
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Figure 7. Unrealistic Body 



Figure 8. Unrealistic Body with Wake 
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Figure 9. Intermediate Body Shape 



Figure 10. Intermediate Body Shape with Wake Body 
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Figure 11. New Body with Minimum Drag 



Figure 12. New Body and Wake with Minimum Drag 
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